Dilepton Enhancement by Thermal Pion Annihilation in the 

CERES Experiment 

H.-J. Schulze and D. Blaschke 
MPG-AG "Theoretische Vielteilchenphysik" , Universitat Rostock, D-18051 Rostock, Germany 



Abstract 



SO 

OS 

On 



We compare the recent CERES data on dielectron production in 200 
GeV/u S+Au collisions with the theoretical predictions due to pion annihi- 
lation in a thermal hadronization and a string fragmentation scenario. Both 
\Q ■ models yield similar results for the dilepton mass spectrum. A satisfactory 

description of the experimental spectrum requires an in-medium reduction of 
the rho-mass and a freeze-out temperature of about 150 MeV in the thermal 
model. We emphasize and discuss the influence of experimental acceptance 
' and resolution corrections. 
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The recent results of the CERES collaboration on the production of electron pairs in 
S+Au collisions at 200 GeV/u |1]|| allow for the first time a quantitative comparison with 
theoretical models for this type of ultrarelativistic heavy ion collisions. The CERES collab- 
oration reports a five-fold enhancement of the dilepton spectrum in the mass range from the 
two- pion threshold up to 1.5 GeV, compared to the expectation based on an independent su- 
perposition of nucleon-nucleon collisions. Their interpretation of this enhancement is being 
due to the production of dileptons from pion-pion annihilation in the hot and dense medium 
that is supposedly formed during the collision. The shape of the observed spectrum then 



hints to a modification of the pion annihilation form factor in this medium. There have been 
a couple of quantitative analyses of the results based on relativistic transport models |3] |5| . 
These works are able to describe the data reasonably well using a suitable medium-modified 
form factor. 

On the other hand, thermal models have been used successfully to describe observables 
like hadron abundances [fj and hadron spectra @ in ultrarelativistic heavy ion collisions. 
While giving relatively good fits to the data, they have the advantage of simplicity and 
and a small number of adjustable parameters, the most prominent one being a "freeze-out" 
temperature Tf. 

In this paper we would like to analyze the dilepton production within a similarly sim- 
ple thermal hadronization model, the standard "benchmark" scenario based on the Bjorken 
model of a boost-invariant longitudinal expansion ||, including a quark phase (u+d+g), a 
mixed phase, and a hadronic pion gas phase. It is characterized by three temperatures, the 
formation temperature Tj of the quark phase, the critical temperature T c for the hadroniza- 
tion in a mixed phase of quarks and pions, and the freeze-out temperature Tf. This scenario 
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has often been employed to predict the dilepton production rate in ultrarelativistic heavy 
ion collisions f|,10], until recently however without a chance of quantitative experimental 



comparison. 

A first test of the model will be the comparison with a purely hadronic scenario based 
on the pion distribution functions as obtained from a string-fragmentation event generator 
(in this case VENUS) ||11|| . Both models should yield similar results for the (unmodified) 
dilepton mass spectrum, because at present collision energies the quark annihilation contri- 
bution in the thermal model is negligible compared to the pion annihilation part 0. This 
supposition is indeed confirmed by our results. The models therefore provide a reliable basis 
upon which more sophisticated features like medium-dependent properties of pion |T2[ or 



rho- meson |T3HT5|, finite pion |T6Hl8| or quark fl9| chemical potentials, transverse expan 



sion [[[J],^(J, or more exotic possibilities [^TJ, can be studied. We will then use the thermal 
model in order to test the influence of a medium modification of the pion form factor. 
With the simple model of an unmodified form factor in the hadronic phase and a dropping 
rho-mass in the mixed phase we are able to achieve a very satisfactory fit of the dilepton 
spectrum when the thermal history of the hot matter is described by a set of temperatures 
(T h T c , T f ) = (250, 160, 150) MeV. 

Before a detailed description of the models, we start however with a discussion of the 
general "kinematical" conditions of the experiment. Obviously, in any comparison with 
a theoretical model, the same kinematical cuts and detector resolutions as in the actual 
experiment have to be applied. In this case, they are fairly severe: A dilepton is only 
detected, if both the electron and the positron momenta satisfy the kinematical restrictions 
on pseudorapidity, transverse momentum, and opening angle, 

2.1 < r] e < 2.65 , (la) 
pf > 0.2 GeV , (lb) 
6 ee > 35 mrad = 2° . (lc) 

These conditions translate into an acceptance function A(M, q?) depending on the mass M 
and transverse momentum of the virtual photon that decays into the lepton pair. A first 
guess can be made by neglecting the transverse momentum. The acceptance for an isotropic 
decay of the virtual photon in its rest frame (as appropriate for a thermal production model) 
is then simply given by 

A(M, )*9(M- AGe W° SlCO ^ n)] ■ (2) 

The numerical value of the second factor is 0.136, and the observed spectrum is therefore 
reduced by nearly an order of magnitude. The full acceptance function can be computed 
numerically and is shown in Fig. [TJ. (In fact it turns out that for two-body decays the 
restriction on the opening angle Eq. (|lc|) has only an effect for extremely low masses M < 
50 MeV, as displayed in the figure). It can be seen that there is a strong variation with 
transverse momentum, in particular for low masses that we are interested in. The notion 
of an acceptance function A(M) depending only on the dilepton mass M, is therefore a 
strong approximation. For a proper treatment the theoretical spectrum dN ee /d 4 q has to be 
multiplied with the acceptance function A(M, g^), before the mass spectrum dN ee / dMdy 
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can be computed by integration over q?. Before a comparison with the experimental data, 
the theoretical spectrum has still to be folded with a Gaussian representing the finite mass 
resolution 5M/M of the detector, which is approximately given by the formula 0: 

5M_ _ f 0.08 , M < 0.5 GeV 

M ~ \ 0.03 + O.lM/GeV , M > 0.5 GeV . ^ 

Finally, for a careful comparison it is important to note that the CERES result represents an 
average over an ensemble of events with mean charged particle multiplicity (dN^/dy) « 125 



and RMS value yj ((dN ch /dy) 2 ) w 135 [TJ. 

We come now to the detailed description of both scenarios of pion annihilation in order 
to compute the dilepton spectrum dN ee /d A q. Starting with a short review of the thermal 
model, the production rate of electron pairs with four-momentum q from the relativistic 
quark or pion gas at temperature T is 

dN " - M ^ {M \f(M,E,T) (4) 



d 4 xd A q 47r(27r) / 



with M 2 = q 2 , a = Jl — 4m 2 /M 2 , {m = m q ,mn), and the elementary cross sections in the 
quark (Q) and hadronic (H) phases, 




a H (M) = \F„(M)\ 2 1 - 1 + —z U 1 - -£ (5b) 



with the pion form factor 



, , XI , Cm A + m 2 V 2 

F n (M) 2 = " ' p 6 

(M z — m l p ) 1 + m-il l p 

and C = 1.3, m p = 760 MeV, T p = 135 MeV. (These parameters fit the experimental data 



22[ better than the standard choice C = 1, m p = 770 MeV, T p = 160 MeV). 

The function F(M, E,T) in Eq. (H) is given by a phase-space integral over the Fermi 



(7a) 
(7b) 



(Bose) distributions of the reacting quarks (pions) |T0|jri 

COsh /3uV ax + COsh Pfi q 



F Q (M,E,T) = i= 1 log 

pq exp \pE) + 1 



F H (M,E,T) = 1 log 

pq exp (pi? — zpfi^) — 1 



COsh /Sc^min + COsh f3jJb q 

cosh/3(a; max - fj^) 



coshp(u min -//„.)-! 



Here E = q ■ u(x) is the dilepton energy, and g = \J E 2 — M 2 is the three-momentum 
component of q in the production fluid element moving with four-velocity u. ^ m ax,min = 
(E ± aq)/2 is the maximal (minimal) quark or pion energy in this frame. The formula 
includes finite chemical potentials for the quarks |T9[ and the pions JTEUTj] ], although we 



restrict ourselves in the following to the case \i q = fi n = 0. The integration of the rate over 
space-time then gives the momentum spectrum of lepton pairs: 
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where r = 37/3 is the ratio of degrees of freedom in the quark and hadronic phases, and G 
represents the integral over space-time rapidity: 



G{M,M T ,T) 



dr) F(M, E,T) 



E(t])=Mt cosh(»7— y) 



(10) 



It sums the dileptons produced in different slices of space-time rapidity r\ at given proper- 
time r (and temperature T). The energy E of the dilepton in these local production frames 
is related by a Lorentz-boost to the energy in the lab system, Mt cosh?/. The normalization 
of the spectrum can be linked to the square of the pion rapidity density by employing the 
relation 7tR 2 tT 3 = ndN^/dy, with k ~ 0.22 in the quark phase. The final result is 



dN^ oi 2 (ndN^/dyf 
~¥q~ " 



167T 



R 2 



7 x 10" 8 GeV 2 



(11) 



where the expression in curly brackets is the one given in Eq. (|9|). In the last step we have 
inserted the appropriate numerical values R ~ Rs ~ 18 GeV -1 , and (dN 7T /dy)nMS ~ 200 @. 

The spectrum is now completely determined by specifying the three temperatures Tj, T c , 
and Tf, and we use here for the moment a standard set of Tj = 250 MeV, T c = 160 MeV, 
and Tf = 120 MeV, that we will discuss later. The comparison with the experimental mass 
spectrum can then be done by multiplying the theoretical spectrum with the acceptance 
function, integrating over the transverse momentum, and folding with the mass resolution of 
the detector, as discussed above. The result is shown in Fig. |2], where we display the spectrum 
before any correction, after acceptance correction, and the final result. The comparison of 
the three curves demonstrates the importance of the corrections. In particular the slope 
of the mass spectrum beyond the rho-peak seems to be entirely determined by the finite 
mass resolution. Before a further discussion of the figure, we turn however to the second, 
completely hadronic scenario of dilepton production, namely by pion annihilation in a string 
fragmentation model. 

We start by repeating the expression for the dilepton spectrum found in Ref. JTT 



dN 



16a 2 



A 27(p T ) 4 (m T ) I 



o 



00 dN n 

dr ^ L (y,r)e 7r (y,r) 



dy 



cosh (y) 



a H {M) exp(-M T /a) 



;i2) 



with d 4 q = MdMdyd 2 qx- This result was obtained by assuming that throughout the reaction 
the pions are produced with a momentum distribution dN n /dydp^ ~ 5(y — r/) exp(— mx/a). 
The parameters (pt) and (m^) are the average transverse momentum and transverse mass 
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of the pions. For further information on the assumptions of the model, and details of the 



calculation, we refer to Ref. [11 



The dynamical content of Eq. (|T2| ) is the one- dimensional integral of the tempo- 
ral development of the product of particle density dN T /dy and energy density e T = 
(hit) /(mR 2 ) dN^/dy in a given slice of rapidity, 

dr -^—(y,T)e 7 ,(y,r) . (13) 

This integral contains the dependence on energy and mass number of the colliding nuclei. 



In Ref. [JTTJ] we used the results of the VENUS event generator |23j on particle and energy 
densities for central O+Au collisions at 200 GeV/u |24]] and found an upper limit of 40 GeV 3 
for this integral. Unfortunately, new results of VENUS on these observables for the case of 
a S+Au reaction are not available. We therefore employ here a simple scaling assumption 
and multiply the O+Au result with a factor 

(14) 



(dNjdy) ) R 



s 



with (dN n /dy)o ~ 150 and (dN n /dy)s ~ 200 appropriate for the two cases. Considering 
the accuracy of the whole model, this procedure should certainly be adequate. Further we 
use in Eq. (|i~2"D a kinematical factor cosh (y) m 2.1, since the center of the CERES rapidity 
coverage lies about 0.6 units away from the nucleon-nucleon cms rapidity (3.0). The resulting 



dilepton spectrum, with an average transverse momentum of pions (px) ~ 0.4 GeV |[25|| , and 
after acceptance and resolution corrections, is also shown in Fig. 

It is immediately noted that although the two scenarios considered here, thermal and 
string fragmentation, rest on completely different physical assumptions and operate on dif- 
ferent time scales, their predictions, at least for the mass spectrum of dileptons, are nearly 
undistinguishable (apart from the quark phase contribution in the thermal model at very 
low masses, which is however completely negligible in the present type of reaction). This 
shows that on the one hand the dilepton mass spectrum is unfortunately rather insensitive 
to the underlying hadron (pion) dynamics, on the other hand it allows us to use now the 
thermal model with some confidence in order to adjust parameters and test assumptions 
about medium modifications. 

Concerning the influence of the choice of temperatures T i; T c , Tf, we find that a variation 
of Tj has practically no importance, since the quark part of the spectrum is negligible anyway. 
The value of the critical temperature T c = 160 MeV appears generally accepted, and we 
fix it as a "canonical" value. A variation of Tf within a relevant range 100 ... 150 MeV 
changes the spectrum by about a factor two, but preserves its overall shape. A low value 
of Tf = 120 MeV, as used in Fig. 0, is a more "traditional" choice PJTOl, whereas rather 



large values are supported by recent fits to hadron abundances in the present experiment 
(Tf « 160 - 170 MeV, g]) or to slopes of hadron spectra (Tf « 150 MeV, f7j). 

We find then that with Tf = 120 (150) MeV and using the free pion form factor the 
theoretical predictions exceed the CERES data by about a factor four (three) in the rho- 
mass region. In comparison, in Refs. M and |J that are based on more sophisticated 
transport models, the corresponding enhancement factors are approximately three and two, 
respectively. We can now conclude that a sizable finite chemical potential of the pions 
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|16|jr7|l can probably be excluded by our result, since this would again enhance the dilepton 



spectrum by nearly an order of magnitude, as shown in Ref. [T? . 

A possibility of reducing the theoretical expectation is of course the in-medium modi- 
fication of the elementary cross sections and form factors, in particular for the rho meson 



26| , |2"T|| . The pion is thought to be rather stable against medium effects ||12|| , and we do not 



consider a modification here. Unfortunately, the numerous models on the modification of 
the rho meson JT73| do not even agree on the simplest qualitative aspects of the variation 
of, say, the rho-mass m p and width T p in the medium. 

We do not intend to introduce a further model here. Instead we will use two simple 
parametrizations of the temperature dependences of m p and T p in the form factor Eq. (0) 
together with the thermal scenario described above in order to assess the principle effects 
on the dilepton spectrum. 



The first parametrization is the one employed in Ref. [pq| , namely an increase of the 
width with temperature, while the mass is kept constant: T p (T) = T p /(1 — T 2 /T 2 ). This 
means that in the mixed phase the pion-pion cross section is reduced to the Born term, while 
in the hadronic phase it is still substantially suppressed. This is probably unrealistic, and 
the results should therefore be considered as a lower bound on what a medium modification 
of the form factor can achieve. 

The second parametrization is based on the assumption that the vector meson mass does 
not change dramatically in a hot hadronic medium and is during the hadronic phase given 



by its vacuum value ||14j| . However, in the mixed phase at the critical temperature for the 



deconfinement transition T c = 160 MeV, the mass might be lowered as conjectured by a 
number of authors ||15|| . Here we simply assume a mass of m p (T c ) = m p /2 « 380 MeV. 
Although at present there is no firm theoretical foundation for such an assumption, we 
consider it as a guess which allows to reproduce the two-peak structure of the experimental 
spectrum. 

We display the results in Figs. |3|,|]. For a better comparison with the data, we have 
added here the part of the dilepton spectrum due to hadron decays, as given in Ref. [|TJ, 
to our pion annihilation results. Indeed the spectrum below 200 MeV is well explained by 
the Dalitz decays of 7Tq, rj, and u mesons. We observe then that the varying- width scenario 
that we consider as a lower bound, is in good agreement with the overall normalization of 
the data using Tf = 120 MeV, whereas the dropping-mass model at Tf = 150 MeV gives 
a rather good fit to the shape of the spectrum, in fact of the same quality as obtained 
with sophisticated transport models and more realistic medium modifications [§-[§]. This 
shows that the larger freeze-out temperature of 150 MeV is better suited to fit the data, 
in agreement with the fairly large freeze-out temperatures deduced from the fit of hadron 
abundances || in the same experiment. A suitable in-medium modification of the pion form 
factor is thus indeed able to account for the observed spectrum. 

This statement is valid with the exception of the data point at M = 0.25 MeV, i.e., below 
the two-pion threshold, that neither in our nor in previous models [|3|-f| can be reproduced 
in a satisfactory manner. This hints to the possibility of a non-negligible modification of 
the pion dispersion relation in the medium, but also the Dalitz decay of the i] meson is still 
important in this region and should be reconsidered. 

We conclude that pion annihilation can indeed account for the observed excess of dilep- 
tons, as already advocated by the CERES collaboration M. We have presented a detailed 
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adaptation of both a simple Bjorken scenario for the thermal dilepton production, as well 
as of a purely hadronic scenario due to string fragmentation, to the experimental situation. 
These models have the advantage of transparency and simplicity, but can not be expected 
to give perfect fits to the data. Nevertheless, we find that both models give surprisingly 
similar results (which means that the possible existence of a quark-gluon plasma can nei- 
ther be affirmed nor rejected from the present data), and can account for the experimental 
spectrum by assuming a dropping rho-mass in the medium. The justification of such an as- 
sumption remains a theoretical challenge. We have here merely provided a guess that allows 
to reproduce the data. A freeze-out temperature of about 150 MeV is supported by the 
experimental results. The correct implementation of acceptance and resolution corrections 
has been seen to be extremely important. 

In particular with a large freeze-out temperature of about 150 MeV nearly all dileptons 
then stem from pion annihilation in the rather long-lived mixed phase (of duration (r— l)r c ~ 
24 fm/c in the one- dimensional expansion scenario, and with the parameters chosen here). 
This puts forward the question of physical reality of such a phase. First principle, i.e., 
lattice gauge calculations do not yet provide an unanimous answer to this problem ||28|| . It 
is possible that a first order transition does not occur at all, or that the release of latent 
heat is fairly small. In any case there exists a relatively narrow transition region between 
the phases, in which hadronic correlations persist. The simple model of a first order phase 
transition between ideal quark and hadron gases is certainly an oversimplification, and should 
be considered a phenomenological approach to mimic the features of this real situation. The 
same remark applies of course to other features of the computation, such as the treatment 
of collision dynamics or medium modification of hadron parameters, as discussed before. 

We would finally like to add a speculative estimate of the expected dilepton spectrum 
in central Pb+Au (or Pb+Pb) collisions |2"§|] . In both models presented here, the dom- 
inant source of change is a scaling factor (dN n / dy) 2 j R 2 , cf. Eq. (|i~4"P. For example, with 
(dN^/dy)^ ~ 700 we would obtain a rescaling of the curves presented in this paper by about 
a factor 12/3.4 « 3.5. Of course, this scaling factor applies only with identical acceptance 
cuts as in the present experiment, so that in practice a completely new analysis suited to 
the experimental conditions will have to be done. 

We would like to thank A. Drees and Th. Ullrich for discussions and detailed information 
on the setup of the CERES experiment. We thank G. Ropke for a critical reading of the 
manuscript. D.B. acknowledges enlightening discussions with K. Redlich. 
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FIGURES 




FIG. 1. The acceptance function for a virtual photon with mass M and transverse momentum 
qr decaying into an electron pair subject to the kinematical restrictions Eq. (|l|). 
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( + iy 'Au , 200 GeV/u ) 
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FIG. 2. The dilepton mass spectrum. The full lines represent the theoretical predictions in 
the thermal (T) and the string model (S) scenario after acceptance and resolution corrections. For 
the thermal scenario, we plot also the results before resolution correction (dashed line) and before 
acceptance correction (dotted line). The data are taken from Ref. [BJ. (Statistical and systematical 
errors added quadratically) . The curve denoted (H) shows the total hadron decay contribution as 
given in that reference. 
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FIG. 3. The dilepton mass spectrum with an in-medium modification of the pion form factor 
and Tf = 120 MeV. The dashed curve shows the result with a temperature dependent width, and 
the dotted curve with a temperature dependent rho mass, as explained in the text. For comparison 
the full curve based on the vacuum pion form factor is also shown again. The hadron decay 
contribution (thin line) has been added to our pion annihilation results. 
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FIG. 4. Same as Fig. but with a temperature Tf = 150 MeV. 
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